Competing Populations in Flows with Chaotic Mixing 
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We investigate the effects of spatial heterogeneity on the coexistence of competing species in the 
case when the heterogeneity is dynamically generated by environmental flows with chaotic mixing 
properties. We show that one of the effects of chaotic advection on the passively advected species 
(such as phytoplankton, or self-replicating macromolecules) is the possibility of coexistence of more 
species than that limited by the number of niches they occupy. We derive a novel set of dynamical 
equations for competing populations. 
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INTRODUCTION 

One of the classical problems of ecology is the identifi- 
cation of the mechanisms responsible for the coexistence 



istence of species in different ways. 

Thus the original problem changed into finding the 
most relevant mecha nisms which maintain divers i ty in 



particular situations (|Connell, 1978 



of competing species. It is an observational fact that son, 1990|; [Tilman fc Pacala, 1993|; [Bartha et al, 1997D 



Huston, 1979|; |wil 



in Nature numerous species are able to coexist, all com- 
peting for a limited number of resources. This observed 
coexistence is at odds with many classical theories and 
empirical studies predicting competitive exclusion of all 
but the most perfectly adapt e d species for e ach limiting 



factor (Cause fc Witt, 1935; Hardin, 1960). However 



one of the key ingredients in these classical studies was 
the assumption of a homogeneous, well mixed and non- 
structured environment which leads to an equilibrium 
state in the system. Thus, if coexistence is to persist 
over longer time periods, it must have a nonequilibrium 
character. 

The coexistence problem is best illustrated in the case 
of phytoplankton communities as was originally presented 
by Hutchinson (1961). Here a number of species coexist 



in a relatively isotropic or unstructured environment, all 
competing for the same sorts of materials, and the num- 
ber of species exceeds considerably the number of limit- 
ing factors. To solve this so called "paradox of plankton" , 
Hutchinson put forward the idea that seasonal environ- 
mental changes prevent competitive exclusion in natu- 
ral phytoplankton communities. Thus the species of the 
community, at least on the time scale of ecological obser- 
vation, are in nonequilibrium coexistence. 

Since then numerous investigations revealed many dif- 
ferent mechanisms, including spatial and temporal het- 



Despite the vivid debate in this field of ecology, there is 
by now a consensus that climatic periodicities and fluctu- 
ations play the main role in causing species' persistence in 



phytoplankton communities (Gaedeke fc Sommer, 1986 
Reynolds, 1993; Sommer et al., 1992 It is frequently 



argued that an intermediate disturbance (Connell, 197?) 
is the most adequate hypothesis for the explanation of 



high diversity in aquatic systems (cf. Reynolds, 1998) 



One can meet a similar problem in early evolution of 
life. Since life evolves from the simple structured enti- 
ties to the most complex ones, there must have been a 
stage in the evolution, when life was essentially no more 
complex than wh at a collection of self-replicating nu - 
cleic acids present ( Maynard Smith fc Szathmary , 1995 ) . 
They were competing for a few limiting resources (such 
as mononucleids and energy rich chemicals) and making 
copies of themselves without any specific enzyme. With- 
out enzymes the copying accuracy could not be very high. 
Estimating the selective superiority of the best replicator 
and the copying accuracy per nucleotide, it is concluded 
that the maxim um length o f these molecules is about 
100 nucleotides (Eigen, 1971). However, in a well mixed 
homogeneous environment, as the prebiotic ocean is sup- 
posed to have been, there are only a few winners of the 
selection, namely the most fit macromolecule surrounded 



erogeneity of habitat, predation, disturbance, coevolu- 1979). But how can we surmount the gap between these 



tion, etc. ( [Wilson, 1990^ [Chesson, 2000] ), increasing the 
probability of competitive coexistence. Naturally, under 
the word "competition" many different biological phe- 
nomena are collected together, which influence the coex- 



by it s closest mutants (Eigen, 1971; Eigen & Schuster 



primitive replicators with 100 nucleotides and the most 
simple RNA viruses with 4000-5000 nucleotides? Speciflc 
replicase enzymes are needed to increase the copying fi- 
delity, and thus the length of the replicator, but these 
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replicators are too short to code specific enzymes 
is the 
Smi 



This 



'Catc h 22" of the prebiotic evolution ( Maynard 
1983): no genome without an enzyme, however 



no enzyme without genomes. This problem can be re- 
solved if some mechanism maintains the coexistence of 
several different replicator molecules and therefore the 
information necessary for coding a replicase enzyme can 
be stored by the union of smaller information carriers. 
In this situation the replication error does not grow ex- 
ponentially as in the case of a base-by-base copying, it 
grows only linearly with the number of smaller carriers. 




FIG. 1. SEAWIFS image of a phytoplankton bloom at 
Shetland Islands, May 12, 2000, from the NASA archive. 
Plankton individuals (light grey) move along a fractal set. 



Current theories p oint out coexistence of replicators 
moving on a surface (Boerlijst & Hogeweg, 1991; Czaran 
& S:;athmary, 2000), preferring thus the concept of "pre- 
biotic pizza" against the concept of "prebiotic soup" 
( Wachtershauser, 1994 ). However some kind of coopera- 



tion among the replicator molecules are built into these 
models, consequently they are not completely competi- 
tive. An alternative explanation assumes that both the 
replicative and enzymatic functions were co-evolved, thus 
the lenght of the replicators and the accuracy of enzy- 
matic functions increased together ( Poole et ai, 1999| 
Scheuring, 2000[ ). In both problems (i.e., in the paradox 
of plankton and in the Catch 22 of prebiotic evolution) 
the traditional population dynamical equation for two 



dNi 

IT 

dN2 
dt 



= a2N2 - 62N2. 



(1) 
(2) 



Here Ni is the instantaneous number of individuals of 
species Bi in a given range of a well strirred region. The 
instantaneous parameters ai , 6i are positive and depend, 
in general, on the concentration of the resource material 
A, too. 
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FIG. 2. The distribution of two populations (light gray (B) 
and dark gray (C)) competing for the same resource material 
(white (A)) in the wake of a cylinder. The flow is from left to 
right. The inset in (a) shows the time-dependence of the pop- 
ulation numbers ns, no and clearly indicates the approach to 
a steady state of coexistence after about 40 time units which 
is the period of the flow. A blowup of the region indicated by 
a rectangle in (a) is seen in (b). Species distribution is strik- 
ingly similar to many of the patterns found on the NASA 
SEA WIFS satelite picture s of plankton blooms (see Fig |l|). 
After [Karolyi et a/. (2000[). 



species B2 competing for the resource A read: 



Independently of the particular form of this depen- 
dence and the dynamical equation of A, no fixed points 
can exist in the system in which both species would be 
in a steady state w ith nonzero values of Ni = N* ( Gur- 
ncy fc Nisbet, 1998| ). The coexistence might, however, be 
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possible in imperfectly stirred environments. As numer- 
ous remote sensing images demonstrate phytoplankton 
are distributed along fractal filaments in the oceans in- 
dicating strong but imperfect mixing environment (Fig 

0). 

Recent development in the field of hydrodynamics en- 
couraged us to revisit the population dynamics of com- 
peting species in open aquatic systems. In aquatic sys- 
tems of large extension, on the time scales characteris- 
tic to the life cycle of microorganisms and replicators, 
the hydrodynamical flows are locally open, i.e. there is 
a net current, transporting both competitors and nutri- 
ents, flowing through the typical observation region. It 
is even more obvious that the flow is open in t he wake of 
islands surrounded by strong ocean currents (Ari'stegui 
et af., 1997) and around the deep see hot spri ngs where 



the cradle of life probably swung (Holm, 1992) 

With the aid of numerical simulations we have pre- 
viously shown that the coexistence of passively ad- 
vected competing species is typical in open chaotic flows 
(Karolyi et ai, 200C; Scheming et al, 200C| ). For simplic- 
ity, we have considered the two-dimensional flow around 
a cylindrical obstacle placed into a uniform background 
flow. For moderate inflow velocities there is a periodic 
detachment of vortices in the wake of the obstacle with 
period T, which forms the von Karman vortex street 
( phariff et al, 19"91 ; lung fc Ziemniak, 1992| ; Sommerei 
et a l, 1996| ^ The flow in the wake is time-dependent 
but still spatially regular. Here individuals of two pas- 
sively advected species compete for a common limiting 
resource, see Fig. ||. We argued that coexistence is due 
to the fractal structures typically appearing in the advec- 
tion patterns of such flows, however, we have given only 
a heuristic interpretation for the mechanisms maintain- 
ing coexistence in this hydrodynamical system. In this 
article we present a mathematical deduction to explain 
coexistence of competitors in open chaotic flows. The 
mathematical problem is to investigate two coupled pop- 
ulation dynamical processes evolving on a fractal sup- 
port. We shall present a new class of equations which 
describe this situation, and allow for the coexistence of 
at least two species competing for the same resource. A 
novel feature of these equations will be a singularly non- 
linear (power-law) form of both the replication term and 
the coupling between the populations. 

In the following section we summarize the qualitative 
features of the relevant physical process, followed by a 
study of the dynamics of a single population in an open 
chaotic flow. Consequently, the coexistence of competi- 
tors is discussed by first giving a qualitative argument 
based on the single population picture, followed by a de- 
tailed mathematical model leading to the aforementioned 
new type of population dynamical equations. Next, this 
theory is compared with further numerical results carried 
out on a simple map modelling the advection dynamics, 
on the so-called baker map. We conclude with a sum- 
mary and outlook. 



PASSIVE ADVECTION IN OPEN FLOWS 



Chaotic advection in open hydrodynamical flows is an 
ubiquitous phenomenon. A flow is considered locally 
open if there is a net current flowing through the ob- 
servation region (Lamb, 1932). It became clear in the 



last decade that passive advection e ven in simple time- 
dependent flows is typically chaotic (Pentck et al.. 1996 
Sommerer et al., 199(\ ; Karolyi & Tel, 1997) and possesses 



complicated particle trajectories. These flows, character- 
ized by strong imperfect mixing, lead to a fractal spatial 
distribution of advected particles in a finite region of the 
flow. This region is called the mixing region. In our 
terminology, a flow is chaotic if the advection dynamics 
generated by the flow is chaotic. 

In the case of several (three or more) types of pas- 
sively advected tracers, di stinguished for examp le by 
their color, it was shown ( Toroczkai et al., 1997 , and 
references therein) that their distribution may follow a 
rather non-trivial topology on the fractal, a property 
called the Wada property: every point on the fractal is 
lying on the boundary of at least three colors. 

It is worth emphasizing that a complicated flow field 
(turbulence) inside the mixing region is not required for 
the flow to be chaotic, (i.e., for complex advection dy- 
namics or for the appearance of fractal patterns). Even 
simple forms of time dependence, e.g. a periodic repeti- 
t ion of the v elocity field with some period T, is sufficient 
( Aref, 1994 ). Thus, for sake of simplicity, we examine 
advection in time-periodic open flows. 

The complicated form of trajectories implies a long 
time spent in the mixing region. In other words, ad- 
vected particles can be temporarily trapped there. It is 
even more surprising, however, that there is an infinity 
of special noncscaping orbits. The simplest among these 
orbits are the periodic ones with periods that are inte- 
ger multiples of the flow's period T. All the nonescaping 
orbits are highly unstable, of saddle type, and possess a 
strictly positive local Lyapunov exponent (which is the 
expanding eigenvalue of the unstable periodic orbit). An- 
other important feature of these orbits is that despite 
their infinite number they are rather exceptional so that 
they cannot fill a finite portion of the phase space. In- 
deed, the union of all nonescaping orbits forms a fractal 
"cloud" of points on any snapshot. This fractal cloud 
moves periodically with the flow and never leaves the 
mixing region. 

Typical advected particle trajectories are not in the 
set of the nonescaping orbits, but are, nevertheless, in- 
fluenced by them. They follow closely some periodic orbit 
for a while and later turn to follow others. This wander- 
ing amongst periodic (or, more generally, nonescaping) 
orbits results in the chaotic motion of passively advected 
particles. Indeed, as long as the particles are in the mix- 
ing region, their trajectories possess a positive average 
Lyapunov exponent A. Hence the union of all nonescap- 
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ing orbits is called the chaotic saddle. The flows relevant 
from our point of view can be conisdered to be incom- 
pressible. This results in a time-reversal invariant, area 
preserving particle dynamics. Therefore, the negative av- 
erage Lyapunov exponent is exactly —A, and it charac- 
terizes the compression towards the chaotic saddle. 

While many of the particles spend a long time in the 
mixing region, the overwhelming majority of them leaves 
this region sooner or later. The decay of their number 
in a fixed frame is typically exponential with a positive 
exponent k (< A), which is independent of the frame, 
i.e., N{t) = A^(0) exp(— This quantity k is the es- 
cape rate from the saddle (or from the mixing region). 
The reciprocal of the escape rate can be considered as 
the average lifetime of chaos, and therefore the chaotic 
advection of pas sive partic les in open flows is a kind of 
transient chaos ( Tel, 1990| ). 

The chaotic saddle is the set of nonescaping orbits 
which advected particles may follow for an arbitrarily 
long time. Each orbit of the set, and therefore the set as 
a whole, has an inflow and an outflow curve, also called 
in the mathematical jargon of chaos theory the stable 
and unstable manifolds, respectively. The inflow curve 
is a set of points along which the saddle can be reached 
after an infinitely long time. The outflow curve is the set 
along which particles lying infinitcsimally close to the 



saddle will eventually leave it m the course of tmie. By 
lookmg at ditterent snapsJiots ot tlicse curves we can ob- 
serve that they move periodically with the period T of 
the flow. Their fractal dimension Dg (1 < Dq < 2 in 
two-dimensional flows) is, however, independent of the 
snapshot. (The inflow and outflow curves have identical 
fractal dimension due to the advection dynamics' time 
reversal invariance.) 

There is a unique relation between the fractal geometry 
and the advection dynamics , expressed by the relation^ 
(iKantz fc Grassberger, 1985| ; |Hsu et ai, 198^ ; [Tel, 1990|) : 



(3) 



It says that the deviation of the dimension from that of 
the plane is given by the ratio of two quantities character- 
izing the global and the local instability of the dynamics. 
Relation (||) shows that out of the three basic character- 
istics (k, a and Dq) only two are independent. When 
speaking about population numbers in what follows, we 



shall use the escape rate and the fractal dimension as 
independent parameters. In the dynamics of the stripe 
widths (see next section), however, only the average Lya- 
punov exponent appears. 

The outflow curve plays a special role since it is the 
only set which can be directly observed in an experi- 
ment. Let us consider a droplet (ensemble) of a large 
number of particles which initially overlaps with the in- 
flow curve. As the droplet is advected into the mixing 
region its shape is strongly deformed, but the ensemble 
comes closer and closer to the chaotic saddle as time goes 
on. Since, however, only a small portion of particles can 
fall very close to the inflow curve, the majority does not 
reach the saddle and starts flowing away from it along its 
outflow corve. Therefore, in open flows droplets of par- 
ticles trace out the outflow curve of the chaotic saddle 
after a sufficiently long time of observation (in fact, the 
populations in Fig. ^ are distributed along the outflow 
curve of the chaotic saddle present in the wake). 



DYNAMICS OF A SINGLE POPULATION 

In this section the mathematical derivation of the dy- 
namics of a single pop ulation living in an open chaotic 
flow is br i efly repeated (Toroczkai et ai, 199^ ; Karolyi ei 
ai, 1999| ; |Tel et ai, 200C| , for more detafls see). Repli- 



cation, competition for the limiting resources, and spon- 
taneous decay are taken into account in our population 
model, while stage and age structure is neglected for sim- 
plicity. We derive discrete and continuous-time models 
as well. 

First we assume that the intake of resource, multiplica- 
tion and decomposition are instantaneous and take place 
at integer multiples of a time lag r. Here r acts as an 
average time-scale on which the reproduction takes place. 

The basic observation is that after a sufficiently long 
time, the filaments of the outflow curve are covered in 
narrow stripes by individuals of species B due to their 
replication (Toroczkai et ai, 1998; Karolyi et ai, 1999). 
Individuals are thus distributed on a fattened-up frac- 
tal set. On linear scales larger than an average width 
e* the distribution of i? is a fractal of the same dimen- 
sion Dq as the outflow curve of the chaotic saddle in 
the flow with passive advection without biological activ- 
ity. Let ef^"^ denote the average width of these stripes 



^By characterizing the dynamics by one single dimension Dq, 
we have assumed that the advection process has a monofrac- 
tal geometry. In reality, a set of dimensions Dq is required 
for the full description of the fractal aspects. It is for the 
q = 1 dimension, the so-called information dimension, Di, 
for which is an exact equality. In practice, however, the 
relative difference between Dq and Di is on the order of a 
few percents and therefore the use of a single dimension is 
justified for practical purposes. 
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right before replication and decomposition takes place. 
It is worth measuring this width in the unit of a char- 
acteristic length scale of the flow (e.g. in the cylinder 
radius in the example of Fig. ||). Thus, e'"-* is a dimen- 
sionless variable. Since material A is available outside of 
these stripes, replication increases the width with some 
constant distance 7, the replication range, while sponta- 
neous decay due to death of individuals decreases it with 
a distance /z. The net effect of the replication and spon- 
taneous decay is then a broadening of the width by an 
amount proportional to the difference a — j — fi, the ef- 
fective replication range. Thus, e(") + ca. Here c 
is a dimensionless number expressing geometrical effects. 
If the fattened-up filaments do not overlap, then replica- 
tion does occur on the both side of stripes and c~ 2. If 
there is overlapping among some of the fattened-up fila- 
ments, such as in the case of a fractal, then c ^ 2. This 
geometrical factor turns out to be slightly time depen- 
dent due to the pulsation of the flow, but for simplicity 
it can be considered to be constant from the po int of view 
of the qualitative behaviour of the population (Toroczkai 
et a 



., 2001) 



In the next period of length r there is no replication 
and decomposition, just contraction towards the out- 
flow curve. The average contraction factor is exp (— Ar), 
where (—A) is the negative average Lyapunov exponent 
of the advection dynamics. Therefore, the width e("+i) 
right before the next replication can be given as 



(n) 



ca]e 



„-At 



(4) 



This is a recursive map for the actual width of the B- 
stripes on snapshots taken with multiples of the time lag 
r. The solution of converges for n — > 00 to the fixed 
point 



(5) 



- 1 



In the time-continuous limit t — > 0, cr ^ 0, but keeping 
a It = Vr constant, one obtains the differential equation: 



de 
'dt 



Ae, 



which has a steady-state solution given by: 



(6) 



(7) 



Knowing the ^-dynamics and that the individuals ac- 
cumulate on a fractal set in the mixing region, the time 
evolution of the number N oi B individuals in that region 
can be calculated. First, note that the area A occupied 
by species B scales a.s A ^ g^-Do .^^ith Do as the frac- 
tal dimension of the outfiow curve for any box size e not 
smaller than the width e of the i?-stripes. We can thus 
choose^ 



^l/(2-Do)^ 



(8) 



If the linear size of the area occupied by a single indi- 
vidual is eo, we have N = Gq^A, and therefore we can 
rewrite @j or so that it represents an equation for the 
individuals in discrete and continuous cases, respectively: 



and 



with 



^ = -kN + q{2 - Do)vrN-'^, 



-2/(2-Z5o) 




q — ce 

Here (||) has been used, and 



/3 



Co 



(10) 



(11) 



(12) 



Here can be interpreted as the net speed of replication. 



appears as a nontrivial exponent.'^ Since the fractal di- 
mension of the outflow curve lies between 1 and 2, ex- 
ponent/3 is positive. For Z?o = 1 the differential equa- 
tion (|lO|) describes a classical surface reaction along a 
line with front velocity Vr in the presence of escape. For 
1 < Dq < 2 it represents a novel form of dynamical 
equations containing also an enhancing biological activity 
term with a negative power of the area occupied by B due 
to the fractality of the outfiow curve. The less B individ- 
uals are present, the more effective the reproduction is, 
because the resolved perimeter is larger. Consequently, 
in a competitive situation the subordinate species has an 
advantage if it becomes rare compared to the dominant 
species. This balancing mechanism can make coexistence 
possible, as shown in the the next sections. 

As one can see from Eqs. (||,^0|), in both the discrete 
and continuum pictures a steady state is reached after 



^In general, (g) also contains a proportionality constant, 
called the Hausdorff volume. Since this only rescales the con- 
stant q in equation (^, for clarity we took the Hausdorff val- 
ume to be unity. 

^For multifractal flows one can show (Tel et al, 2000) that 
exponent /? is that given by ( [l2| ) with Do replaced by the 
information dimension Di. 
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a sufSciently long time if the geometrical factor c (and 
therefore also q) is constant ( Toroczkai et ai, 2001 ). In 
this case, the steady-state number of individuals in the 
mixing region is N* = e^^(e*)^~^'' where e* is given by 
and (0) for the discrete and continuum cases, respec- 
tively. 



A MODEL OF COMPETITION 



and I is the constant inflow of resource A into the region 
of observation. Equations ( |l4|) correspond to the 

general scheme (||) and (||) given in the Introduction by 
identifying /i^ with Si and ^iA (where A is given by the 
right hand side of dlq)) with a^. 

After analyzing (]lj) and (^5|), one can easily see that 
species with lower ratio "fi/ni of replication and death 
rates would be outcompeted, and thus stable coexistence 
is impossible. 



As in the single species case, we consider a simple 
model of replication and competition with passively ad- 
vected point like individuals of type Bi and B2, multi- 
plying themselves instantaneously. The resource material 
A which the different species Bi and B2 compete for is 
uniformly distributed on the surface of the flow. Sim- 
ilarly to the single species case, the parameters 7^ and 
/Xj (« = 1,2) are defined as the increase and decrease of 
the Bi stripe width due to replication and decomposi- 
tion, respectively. Therefore two autocatalytic processes 
A+Bi — >2Bi,Bi — > A and A+B2 — > 2B2, B2 — > A 
represent the replication and competition process in our 
model in an imperfectly mixed environment. Similarly to 
the single species case, the parameters 7^ and /i; (z = 1, 2) 
are defined as the increase and decrease of the Bi stripe 
width due to replication and decomposition, respectively, 
so that the effective replication distances are ai = ji — fii- 

As before, an important feature of the advection dy- 
namics is its deterministic nature. Concerning the pop- 
ulation dynamics, this implies that we work in the limit 
of weak diffusion and assume that the mutual diffusion 
coefficients between any pair of the constituents is small. 

Prior to discussing the consequences of the imperfect 
mixing generated by the chaotic flow to this dynamics, it 
is worth briefly giving the traditional equations governing 
the above defined autocatalytic processes in a well-mixed 
environment. In a fixed region of observation they are: 



dN2 
dt 



-fiANi - niNi, 



^2AN2 - fl2N2, 



(13) 



(14) 



where Ni denotes the number of individuals of species 
Bi, and A is the instantaneous amount of the resource 
material in the same region. Note that the meaning of 
the replication and death rates are slightly different here 
from those in the discrete model (thus e.g., /i^ in eqs. 
( |l3| ) , ( |r^ ) is of dimension frequency, while the same quan- 
tity in the discrete version is a distance). If the dynamics 
of resource is much faster than the dynamics of compet- 
ing species, then the former can be considered to be in 
a quasistationary state: dA/dt = 0. The equation for 
resource A is then 

dA 

— = = /-7iA7Vi -72AiV2, (15) 



COEXISTENCE OF COMPETING SPECIES 

In the following, using a gedankenexperiment, we find 
conditions under which the two species may coexist in the 
imperfectly mixing environment characterized by the ex- 
istence of a chaotic saddle and its (fractal) outflow curve. 
First, let us assume that initially there is only one species, 
for e.g., B2 in the mixing region. Also for the simplicity 
of the writing, we will refer to the continuum version of 
the single species dynamics, namely Eq. (10). After a 
time the number N2 of B2 will be close to a steady state 
value, A'|. We now let a small quantity of species Bi 
invade the mixing region, so small that it cannot change 
the steady state of B2. The question is: under what 
conditions this invasion will lead to a self-sustained Bi 
population, coexisting with B2"! 

When we are letting species Bi invade the mixing re- 
gion such that B2 is already present, we must ensure for 
the coexistence that there is not only Bi — B2 interface 
present but also Bi — A interface, at all times. To show 
that this is indeed possible, we recall the Wada property 
of mixing on the chaotic saddle, mentioned previously. 
As shown by (Toroczkai et ai, 1997), if species Bi is in- 
serted such that the Bi— A interface intersects the inflow 
curve, then this boundary must be present arbitrarily 
close to all points on the outflow curve. Therefore this ini- 
tial condition ensures the existence of the Bi—A interface 
at all times on the saddle and its outflow curve. Thus, 
when there are interactions among the species, then, due 
to the relatively long time they spend on this fractal, and 
due to the largely increased interfaces a non-trivial and 
novel type of behavior can emerge. 

Species Bi becomes thus trapped by the chaotic sad- 
dle, and will be distributed in very narrow filaments along 
the surface of some of the B2 stripes close to the Dq di- 
mensional outflow curve. Since by assumption, the pop- 
ulation number A^i of species Bi is small, N1/N2 ^ 1, 
there is no feedback on B2 and this species remains in 
dynamical equilibrium: dN2/dt « 0. The dynamics of 
species Bi can thus be described by ( p^ written for iVi. 
The important nature of (^0|) is that the first term on the 
right side, (which is responsible for the outflow from the 
fractal set) tends to zero for very small iVi. The second 
term, however, which describes an autocatalytic process 
can be arbitrarily large if iVi is small, due to the nega- 
tive exponent —f3. Thus, if Ni is sufficiently small then 
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dNi/dt > 0. It means that the Bi population always 
increases if it's number is close to zero. In other words, 
the fixed point = is unstable. 

Following the same argument, if Bi and B2 are al- 
ready in dynamical equilibrium on the fractal set then a 
third species can invade this coalition in the same man- 
ner. This coexistence has indeed been demonstrated for 



three species by numerical simulations (Karolyi et al. 
200(|). 

In this gedankenexperiment we assumed that the in- 
vading population is so weak that it is distributed in nar- 
row filaments along the "surface" of the already existing 
population and does not infiuence this population at all. 
After the number of the invading population has started 
to grow, there is an increasing interaction between the 
populations due to the competition for the same resource. 
As a consequence, the originally steady population is not 
in a stationary state and the dynamics should be treated 
in a self consistent manner. This is shown in the next 
section. 



A MATHEMATICAL MODEL FOR THE 
COMPETITION DYNAMICS 

After sufficiently long time, both species Bi and B2 
will be distributed in narrow stripes along the chaotic 
saddle's outflow curve as follows from the passive advec- 
tion dynamics. Due to the replication and decomposi- 
tion, however, the stripes have finite widths (cf. Fig. ||) 
which might depend on time. Let e" denote the dimen- 
sionless average width of the stripes right before an in- 
stantaneous replication takes place. These stripes are de- 
fined by the fact that outside of them there is only back- 
ground material A available. Inside the stripe of width 
e*^"^ there might be several narrow Bi or B2 filaments. 
The background material A is eaten up sooner or later 
in the inside of any stripe, therefore, for the sake of an 
easier presentation, we assume that this is the case and 
only material Bi and B2 are present. Let us denote the 
total widths of all the filaments of a given material within 
an e*^"^ stripe by e^"'' with i = 1,2 corresponding to Bi 
and B2, respectively. The sum of these partial widths 
is of course the total one e^"' -|- ej""* = s''"^- Our aim is 
to build up the dynamics of the partial widths based on 
plausible assumptions, from which the dynamics of the 
different populations follows. 

We assume, that the boundaries are occupied by 
species Bi or B2 with probabilities pi and p2 , respectively. 
In other words, a stripe-boundary picked at random from 
the many filaments of the outfiow curve will have a prob- 
ability Pi to be of type Bi, i — 1,2. If mixing of the two 
species were perfect along the fractal set, these proba- 
bilities would be equal to their relative densities. This 
is not the case, however. The relative position of the 
species in the initial distribution to the inflow curve de- 
termines which individual or patch of individuals will be 



trapped by which orbit of the chaotic saddle. The rest, 
i.e., the untrapped individuals will drift out of the mix- 
ing region. The trapped individuals, however, will stay 
there forever, and follow their specific trapped orbit. In 
the course of time, individuals give birth to others of the 
same species, and patches of individuals are stretched 
along the outflow curve specific to the trapping orbit of 
the chaotic saddle. In either cases, we end up with long 
stripes of the two species lined up along each other in an 
alternating manner, tracing out the outflow curve. Then 
the probability of one species to be on the edge of these 
lines, and thus to be capable of reproduction, depends 
on which trapping orbit will produce the filament of out- 
flow curve being on the edge of the stripe, on the order 
in which the species are lined up across one stripe, and 
on the actual width of the coverage of the fllaments. In 
other words, it is the complex chaotic dynamics which 
makes the introduction of probabilistic concepts — on a 
somewhat phenomenological level — unavoidable. 

The probabilities pi depend on what the distribution of 
the species inside the stripes is. Thus, the simplest pos- 
sible assumption is that the probabilities depend on the 

(n) 

partial widths . Their actual functional form might 
also contain parameters of the flow and of the biological 
activity. 

(n) 

Naturally the probabilities fulfill < Pi < 1 and 
p^"^ = 1 — p^"'' . They might have a general dependence on 
the partial widths e|"\ i = 1,2: p-"-* = p-"-* (^£i"\ £2"'') ■ 
Due to dimensional reasons they can only depend on the 
ratio z'") = e^"''/e2"''- Thus we write 



(n) 

Pi 



pi-' - 1 



(16) 



Here a; is a parameter of the distributions, and incorpo- 
rates the dependence on the replication rates. We also 
made the plausible assumption, that g has no explicit n 
(or time) dependence. A general property of g is that it 
vanishes in the origin g{0) ~ since this expresses the 
obvious fact that if species Bi is missing, then the prob- 
ability to find it in the filaments is zero. Similarly for 
infinitely large values of z it must be unity: g{oo) = 1 
which corresponds to the absence of B2- Also, due to 
the fact that g{z) is a probability distribution, we must 
have < g{z) < 1 for all z > 0. Furthermore, the func- 
tional form must be symmetric by interchanging the role 
of the species. This implies that one must have ^'(O) > 0, 
where the prime denotes derivation with respect to the 
argument. This implies 



P2 



g{l/z("\l/u;) 



(17) 



where the appearance of 1 /w means that an interchange 
of the species index brings the parameter in its reciprocal 
value, as e.g. in the case when uj = oxj 02 (the dependence 
on the ratio of the replication distances follows from di- 
mensional reasons) . The normalization of the probability 
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implies 



g{z,u)+g{l/z,l/u) = 1 



(18) 



This is a functional equation for g. With the above prop- 
erties and boundary conditions we find that a family of 
solutions is given by the form: 



.9(^) 



(19) 



with a and uj as two positive parameters. In the range 
of < a < 1 the smaller population is less probable 
on the boundary but yet with a weight which is weaker 
than linear in the widths. For a = there is no width- 
dependence at all, the probabilities pi are constant. The 
case a — 1 and uj ^ I corresponds to a homogeneous 
mixing within the stripe of width s. For a > 1 a super- 
dominance is described. In the next section we show that 
the form ( p^ of g{z) is indeed in good agreement with 
numerical simulations, and determine values for param- 
eters a and lu. 

The broadening of the average widths is then cfiipj"'' 

(n) 

and ca2P2 due to species Bi and B2 , respectively. Here 
the geometrical factor c and parameter ai = — /ii have 
the same meanings as in the single species problem de- 
fined previously. 

Thus, similar to (||) the partial width of Bi after the 
(n -I- l)st step is 



.(n+l) 



CCTiP. 



(n) 



(20) 



for i — 1,2. Note that in our theory, cpi and cp2 can 
also be interpreted as renormalized geometric factors for 
each species, due to the screening effects at the bound- 
aries of the stripes. As a consequence, the total width of 
the stripes changes at a replication as 



,(n+l) _ 



.(") 



c{aipl 



(n) 



a2P2 



(21) 



For simplicity, the explicit width-dependence ( [161) of the 
probabilities has not been written out. For cti = a2 we 
recover the width dynamics of the single species problem, 
see(|). 

Next we turn to the dynamics of the number of indi- 
viduals. On scales larger than or equal to e^"^ the total 
number of individuals N = Ni + N2 occupied by stripes 
appears to be a fractal of the same dimension Dq as the 
outflow curve. For simplicity of writing we assume that 
both species have the same size Sq (an extension for dif- 
ferent sizes is straightforward). 

Since the relation between the and the number 
of individuals 7V^"^ is the same as in the single species 
model, we can use (^). Thus ( |2l] ) implies a recursion for 
the area right before replication as 



l/(2-D„) 



(") 



(n) 
(T2P2 



2- Do 



(22) 



with q given by (11). Next, we derive the dynamics of 

(n) 

the number of individuals for species i contained 

in the stripes. The number of individuals of species i is 
the portion of the total number Af(") proportional to the 
partial widths: 



(n) 



(23) 



This is due to the fact that there is no fractal scaling 
below e("). Since = [eg7V(")]i/(2--D«), Eq. (|||) leads 
to 



-(") ^ Ar(») 2/(2-Do) r^(„) 



(24) 



As another consequence of (23), the ratio of the partial 
widths is the ratio of the population numbers: 



.(») 
\i 

-(») 



(") 



N. 



in)- 



(25) 



From Eqs. ( P0[ ) and (24) we therefore obtain the dynam- 
ics of the population numbers as 



r(») 



/3 



for i = 1, 2. Here exponent /? is the same expression ( [l^ ) 
as in the case of the single species problem, and q is given 
by (0). 

We observe that by dividing the rearranged ( p6| ) for 
i = 1 by that with i = 2, one obtains 



(n) 

in) ' 



(27) 



where 



~ e^iV^"+'^ [iV("+i)]/3 - 7V^"^ [7V(")]/3 ■ (28) 



From this, pi = 1 — p2 is easily found as 



(") 
Pi 



M(") 



(29) 



This relation provides us with a method for measuring 
how the probability pi{z) = g{z) depends on the ratio 
z = N1/N2 at any instant of time. We shall use this 
observation to extract the form of the g function from 
numerical results. The ratio of fixed points, z* = can 
be calculated from (p9|) assuming that N2 ^ as 



M* 



(72 



(30) 
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The time continuous limit is obtained by letting both 
the time lag and the effective replication ranges go to 
zero so that their ratios remain finite. Thus we define 
replication velocities 

Vi = lim — , (31) 
T^o r 

with 1 = 1,2 for species Bi, B2, respectively. In the con- 
tinuous time limit, the differential equations obtained for 
the partial widths from ( pO| ) read as 

= -\ei + cvipi{ei/e2) (32) 

where pi — g, p2 = ^ — g- 

The differential equation for the number of all individ- 
uals follows from ( p2| ) as 

— = -kN + q{2- Do)vN~^. (33) 

Here 



If one of the species, say i32, is not present, thenpi = 1, 
-^2 = P2 = ^'2 = and hence v = vi = Vr, N = Ni and 
Eq. ( |33| ) becomes equivalent to ( [To| ) . The same happens if 
both species are equivalent, i.e., for vi — V2 when v = Vr- 

A simple further equivalent form can be derived by 
using relative densities = Ni/N . The equations de- 
scribing the populations then become (by using (|3^ ) and 

^ = qN-^-^{vipiC2 - V2P2C1) (37) 

with ci + C2 = 1. The temporal change of the densities 
is determined by the weighted relative difference in the 
densities. Note that they are muUiplicatively coupled to 
N~^~^ which is proportional to the average width of the 
filament covering. For Dg = 1 this is just 1/N and it 
is the spatial concentration or the density of the total 
population. For 2 > Dq > 1 (fractals), this factor is the 
fractal spatial density of the population as a whole. 



: Pivi +P2V2 



(34) 



is an average velocity, but note that it is not a constant 
since the pi depend on the population numbers, and q is 
given by (|ll]). 

The differential equation for the number iV^ of individ- 
uals of the two species can be derived from ( p2| ) and the 
continuum version of (p4[), i.e., — e^^^^ ^""^ NiN^^ . The 
result is 

^ = ~kN, - q{Da - l)vN-^-^N, + 
at 

qv,p,{N^/N2)N-P, (35) 

with TV = iVi + 7V2. Here (||) and (|l2|) have been used. 
By summing over i in ( ^5| ) one recovers Eq. (|3^). 

An equivalent form is obtained after rearranging terms 
and taking into account the definition of the average 
replication velocity (pj). It reads 



IT 



-kNi + qN-f^'^ [(2 - Do)vNi + 
{V1P1N2 - V2P2N1)] , 



(36) 



and an analogous expression for the second species ob- 
tained from (^) by interchanging the indices 1 and 2. It 
can be clearly seen that the first term of the bracket cor- 
responds to the growth of the total population, while the 
second describes the effect due to a weighted difference in 
the population numbers. Expression ( p5| ) or ( ^ ) repre- 
sents a strongly coupled set of nonlinear equations with 
a novel type of power-law behavior (with negative expo- 
nent —f3). This set of equations is the central result of our 
paper since it can be considered as a population dynam- 
ics describing the coupling of two populations mixing on 
a fractal, and as we show below, opens up the possibility 
to have a nontrivial coexistence. 



COEXISTENCE ANALYSIS: FIXED POINTS 
AND THEIR STABILITIES 



It is simple and instructive to study the time- 
continuous dynamics of the widths Si, i — 1,2 in and 
around steady states. Assuming, as before, c = const., 
we find from 



Ae* = cv,p, {el/ 6*2) . (38) 

The weighted sum of the fixed points widths ( ^8| ) gives 

CV1V2 



£lV2 + £2^1 



A 



(39) 



From here on, for the stability analysis, we shall use the 
explicit form (|l9|) for the function g. Thus (^8|) translates 
into: 



^^■2=CV2-^E——^ 



(40) 
(41) 



Formula (^9|) shows that one species always survives. 
Without loss of generality, we can choose this to be 
species B2, and everything remains valid with the indices 
1 and 2 switched. It is worth defining: 



z =eje2. 



(42) 



Since e*2 is not zero, the ratio of the fixed point equations 
(|^,|l|) yields: 



z ~ or z = 

f2W 



(43) 



9 



The first equality describes the z* 7^ coexistence fixed- et ai, 2000), so we deal only with the quantitative verifi- 



point while the second describes the non- coexistence 
fixedpoint. 

The equations (^) of the continuous case written out 
explicitely are as follows: 



dei 
~dt 
de2 
~dt 



-Asi + cvi- 



- we. 



CV2- 



(44) 
(45) 



The linear stability of a fixed point (e* > 0, £3 > 0) 
will be given by the eigenvalues of the stability matrix 
E, calculated from (H]45|) as follows (here we also used 

(to: 



E = A 



CV2 ^ CV2 -L 

oA 1 I aA 

CVi ^ CVi i-/ 



The eigenvalues of E are easily calculated: 



A+ = -A(l - a), A_ = -A. 



(46) 



One eigenvalue of the width dynamics is always the 
negative of the chaotic advection's positive Lyapunov ex- 
ponent. As long as the parameter a is less than unity, 
the other eigenvalue is also negative. 

We find that for < a < 1 coexistence is stable, for 
a > 1 it becomes unstable. 

The case a = 1 is special. It follows from (^3|) that 
for a — 1, and lu ^ vi/v2, the non-coexsitence point 
is the only fixed point, and it is stable. (If uj > vi/v2 
then {el > 0,^2 = 0) is stable, if u < Vi/v2 then 
(e* =0,63 > 0) is stable fixed point.) Having a = 1 
with ui = Vi/v2 implies that all points fulfilling (|39| ) are 
fixed points of marginal stability. Thus, stable coexis- 
tence is found in the 



< a < 1 



(47) 



regime. 

Interestingly, the stability conditions are the same for 
discrete-time dynamics, as well. 

In the next section we analyze coexistence in a sim- 
ple chaotic dynamical system, the Baker map, where we 
show that g{z) is given by Eq. ( pj| ) in this process, and 
determine the parameters a and u from numerical exper- 
iments. 



NUMERICAL RESULTS 

In this section we present numerical verification of the 
new type of population dynamical equation we intro- 
duced before. We have already shown that coexistence 
in open flows is possible (Scheuring et ai, 2000; Karolyi 



cation of the theoretical results. 
For computational simplicity. 



we use the so-called 



baker map to model the flow (Toroczkai et at, 2001). 
This can be considered as a simplified model of stretching 
and folding in a chaotic flow observed periodically after 
specified time-intervals. Thus, in this case r, the time lag 
between instantaneos multiplications of the species, is an 
integer number denoting the number of snapshots taken 
of the flow between two consequtive multiplications. The 
baker map, acting on the unit square, gives the new lo- 
cation {x',y') of an individual starting at point {x,y): 



ax + {l-a)0{y-l/2), xe[0,l], 
iy-(i-l)0(y-l/2), ye [0,1], 



(48) 



where a < 1/2 is the parameter of the baker map, and 
6{x) is the Heaviside step function. The action of the 
baker map is shown schematically in Fig. ^. The area 
preserving property of this baker map models the incom- 
pressibility of realistic hydrodynamical flows, while out- 
flow is modelled by neglecting the area hanging over the 
edge of the unit square. 
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FIG. 3. Two consecutive steps of the baker map and two 
replications (r = 1) for the single species model. The bands 
of width a become occupied by B in each replication. The 
material hanging over the unit square is discarded. 

Starting with any initial conditions, after a few steps of 
iterations both species will be distributed along narrow 
filaments parallel to the y axis. 
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FIG. 4. Reaching the equihbrium states with the coexis- 
tence of two species is shown. Initially, two patches of species 
were placed, one patch of Bi in x £ [0;0.1], y G [0;1], and 
another patch of B2 at x £ [0.1; 1], y £ [0; 1]. The param- 
eter values are a = 0.4 for the baker map, and cti = 0.003 
(72 = 0.001 for the competing species. The areas covered 
by the species is shown right after the multiplications taking 
place. After an initial transient (time-steps 1-4), we have a 
rapid convergence to the fixed point (time-steps 5-18), after 
that we have an equilibrium setting in (time-steps 18-20). 



bution g. Figure ^ shows pi as a function of N1/N2 for 
fixed parameter values, but for various initial conditions. 
There is a single function covering the measured points 
which can well be fitted by the form gi{z) = z°'/{z°' +uj. 
In all cases a < 1 was measured indicating that the coex- 
istence fixed point is stable. Also note that u! « (^1/^2)" 
was found in all experiments, which implies that (Hq) 
holds. 



After T baker steps, individuals of species Bi multiply 
and give birth into a vertical stripe of width ai covered 
by resource A, lying along the borderline of the previ- 
ously occupied region of species Bi parallel to the y axis. 
In the numerical experiments, we used t — 1, that is, the 
species multiplied after each baker-step. Regions which 
are invaded by both species after instantaneous multipli- 
cation are divided between them in a ratio of (Ti/a2- It 
is expected that describes the time-evolution of the 
species, reaching the fixed-point ( |30| ) eventually. Figure^ 
shows in a typical case how the equilibrium state with 
coexistence is reached after about 18 baker steps. Sim- 
ilar results were obtained with various other parameter 
settings, in accordance with the theoretical results. We 
checked the validity of the form (^0|) against the numeri- 
cal results in steady states. N* is the fixed-point number 
of individuals of species Bi. The fixed point values are 
found to fulfill (|^, see Fig. |. 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 








0.1 



10 



=Ni/N2 



FIG. 5. The dependence of the probabilities pi on A''i /N2 in 
the nontrivial fixed points. The initial positions do not influ- 
ence the fixed point reached. The curve g{z) = z j (z + a\ j 02) 
is shown with solid line for a\ — 0.002, (T2 = 0.001, with 
dashed line for ai — 0.003, 0-2 — 0.001, and with dotted line 
for (Ti = 0.004, (72 — 0.001. All the measured fixed point 
values fulfill g{z*) = z* /{z* + axja-i). The fixed points are 
marked by crosses (a = 0.25), black squares (a — 0.3), starts 
(a — 0.35), and circles (a — 0.4 as the baker parameter). 

Next we check the validity of (|2^) for the time- 
evolution before reaching the convergence. We measure 
the population numbers in discrete time n and use re- 
lation (p9) to extract the form of the probability distri- 




FIG. 6. The dependence of the probabilities pi on N1/N2 
with (Ji = 0.003, (T2 = 0.001 values, and with five differ- 
ent initial conditions. The parameter of the baker map was 
a = 0.4. The initial positions do not influence the fixed point 
reached. Solid line shows the function g{z) — z" / {z°' +u}) with 
a = 0.818 ± 0.002 and uj = 2.312 ± 0.006. The parameter of 
the baker map was a = 0.4. 



We also measured how a depends on the parameter 
of the baker map, or, on the fractal dimension Dq of 
the outflow curve of the chaotic saddle. We found that 
a = 0.79 Ina + 1.54, see Fig. 0. Using the fact that 
Do = ln2/;77,(l/a), we obtain a = 1.54 - O.SS/Dq. 
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FIG. 7. Dependence of a on the parameter a of the baker 
map. Various points are measured values for multiple ai/a2 
ratios. The dashed line is the function a{a) — 0.79 In a -I- 1.54. 
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DISCUSSION 

We derived a novel type of coupled population dynamic 
equations for two populations competing on a fractal set 
provided by open chaotic flows. The equation for the 
number of species in a given fixed range of the flow can 
be written in the general scheme (cf. |36| ) 



dNi 
dt 



dN2 
dt 



ai 



N ' 1 



(49) 



(50) 



The coefficients ai of the replication terms depend on the 
relative densitites (denoted by Ni/N = Ci) only. Their 
explicit form follows from (pq). For example, 



El 

N 



= '^1 AT 



(l-^o).f 



VlPl 



(51) 



The structure of these equations is similar to that of 
(^, (H) or (|l3|, |l^. The time derivative is the sum of 
a gain term and a loss term, but now the gain term 
contains a nontrivial negative power of the population 
number and is coupled to the other population in a non- 
linear way. These equations describe the population dy- 
namics in an imperfectly mixed enviroment of dimension 
1 < Do < 2. The fractality Dq of the mixing region 
(in our case of the outflow curve) appears in the power 
/3 = (2 — Dq)/{Do — 1). In this set of equations a phe- 
nomenological function {pi) is also present characterizing 
the probability that a given population is on the sur- 
face of the fractal support with free access to the single 
available resource. Based on general arguments and a 
simple model, this function turned out to be a normal- 



ized power law distribution of the type of (19). This 
form expresses a kind of "advantage of rarity" principle: 
for exponent < a < 1 the derivative is infinite in the 
origin, a very small increase in the size of the weaker 
population leads to a drastical increase of the probabil- 
ity for being on the free surface and hence to grow. On 
the contrary, for a > 1, only a relatively large popula- 
tion size has considerable growing probability, in this case 



the weaker population dies out. It is worth mentioning 



that in the gedankenexperiment of Section "Coexistence 
of competing species", we did not see this effect since 
we had an initial stage at which the weaker population 
does not yet influence the stronger one. The probabil- 
ity of being on the surface was assumed to be constant. 
It is the interaction between the two populations which 
lead s to the power law distribution. Its exponent is de- 



tern lined by the flow and the biological process. With 
exponents larger than unity this form does not allow for 
coexistence. The presented mathematical forms and the 
conditions for coexistence remain valid if m > 2 species 
live in open chaotic flow, a numerical evidence for which 



has been reported by Karolyi et al. (20"00| ). It is natural 
to expect that the probabilities pi, i — 1, m appear in 
the generalized form of pi = (t^isf )/(X]i=i ^i^f )' where 
Ei are the partial width of the species and uji are phe- 
nomenological constants. 

Although in the numerical simulation, based on the 
baker map as a model flow, we only found stable co- 
existence, we also carried out simulations where the bi- 
ological process was not based on parallel stripes filled 
out homogeneously with individuals, as assumed in our 
theory. In these cellular automaton-like simulations the 
replication and competition process is carried out on a 
uniform rectangular grid of lattice size eq. This eo can be 
considered as the smallest distance between the individ- 
uals, or the linear size of a single individual below which 
there is hard-core exclusion among them. Individuals 
of each species can occupy the center of each grid-cell. 
When they are adverted by the flow into another grid- 
cell during the time t, they are instantaneously placed 
to the center of that grid-cell. During reproduction, they 
give birth to new individuals in the surrounding empty 
grid-cells, whose centers are within a distance ai. If more 
than one species tries to give birth into the same grid- 
cell, only one of them will be able to do so according to 
one of the following rules: rule I: both species can win 
this competition in each cell with equal probability, or 
rule II: both species can win this competition in each 
cell with probability proportional to the number of indi- 
viduals of the same species intending to give birth there, 
or rule III: always the better competitor (with higher a) 
wins. Our results show that the coexistence depends on 
which rule has been applied. In some cases one of the 
populations was competed out, but even in such cases 
the distribution function was found to be of the shown 
form, with an exponent a > 1, in full harmony with 
the theory. It is worth mentioning that wit h the same 
rules on the lattice, in previous simulations ( Karolyi ei 



ai, 2000) for the more realistic fluid dynamical case of 



a flow around an obstacle we always found coexistence. 
This indicates that the boundary layer present around 
the obstacle enhances the chances of survival. 

Our theory does not describe the effects of diffusion. 
Besides the fact that for individuals of small but macro- 
scopic mass and size, like eg. phytoplankton, diffusion is 



not believed to be important, it can be shown (Tel et al. 



200C ) that weak diffusion in such models only renormal- 
izes the replication rates. As a consequence, the cut-off 
scale below which fractality cannot be observed is some- 
what increased, but the population dynamical equations 
remain unchanged. 

In this theory the location dependence of the death 
and replication rates is not taken into account. Such ef- 
fects call be studied in numerica l simulations ( pK!arolyi~ 



al, 2000| Bantoboni et al, 2001 ) and are not expected to 
change the essence of our findings. 

In conclusion, we have shown that a particle-like (mi- 
croscopic) model of individuals competing for a single 
resource around a fractal outflow curve of a chaotic flow 
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leads, on the level of the total number of individuals, to 
dynamical equations with unusual singular terms. These 
describe enhanced competition due to inhomogeneous 
mixing and a kind of advantage of rarity property. The 
appearance of unusual population dynamical equations 
can be expected in general in all cases where the individ- 
ual dynamics is not taking place on full compact regions 
of the space but are restricted to fractal subsets of it. 
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